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43 ABSTRACT 

Pk 

Aims. To test the dust torus model for active galactic nuclei directly, we study the extent and morphology of the nuclear dust distribu- 
tion in the Circinus galaxy using high resolution interferometric observations in the mid-infrared. 

Methods. Observations were obtained with the MIDI instrument at the Very Large Telescope Interferometer. The 21 visibility points 
recorded are dispersed with a spectral resolution of A/ 6 A ss 30 in the wavelength range from 8 to 13 /im. To interpret the data we used 
a stepwise approach of modelling with increasing complexity. The final model consists of two black body Gaussian distributions with 
dust extinction. 

Results. We find that the dust distribution in the nucleus of Circinus can be explained by two components, a dense and warm disk-like 
^ ' component of 0.4 pc size and a slightly cooler, geometrically thick torus component with a size of 2.0 pc. The disk component is 

^\ ' oriented perpendicular to the ionisation cone and outflow and seems to show the silicate feature at lOyum in emission. It coincides 

with a nuclear maser disk in orientation and size. From the energy needed to heat the dust, we infer a luminosity of the accretion disk 
of Lace = 10'° L , which corresponds to 20 % of the Eddington luminosity of the nuclear black hole. We find that the interferometric 
data are inconsistent with a simple, smooth and axisymmetric dust emission. The irregular behaviour of the visibilities and the shallow 
decrease of the dust temperature with radius provide strong evidence for a clumpy or filamentary dust structure. We see no evidence 
for dust reprocessing, as the silicate absorption profile is consistent with that of standard galactic dust. We argue that the collimation 
of the ionising radiation must originate in the geometrically thick torus component. 

Conclusions. Based on a great leap forward in the quality and quantity of interferometric data, our findings confirm the presence of a 
geometrically thick, torus-like dust distribution in the nucleus of Circinus, as required in unified schemes of Seyfert galaxies. Several 
aspects of our data require that this torus is irregular, or "clumpy". 

Key words, galaxies: active, galaxies: nuclei, galaxies: Seyfert, galaxies: individual: Circinus, radiation mechanisms: thermal, tech- 
niques: interferometric 



1. Introduction allegedly is responsible for the anisotropic absorption and ob- 
scuration has been missing. It is too small to be significantly re- 

In the widely accepted model for low and intermediate luminos- solyed wkh single difjh telescopes in the infraredi even for 10 m 

ity active galactic nuclei (AGN), the central engine, consisting of dass telescopes at the diffraction hmit , It has hen ce been impos- 

a hot accretion disk around a supermassive black hole, and the sible tQ directly confirm or disprove the existence and geometry 

broad line region (BLR) are surrounded by a torus of obscuring of tWs essential compon ent. Properties of the dust distribution, 

dust. Where the radiation from the central engine can escape, an such afj itg ^ temperature> geometry, alignment, and filling 

ionisation cone and the narrow line region (NLR) are formed. factor haye hence remained unknown . A new approach is possi- 

Although motivated by the observation of broad emission lines We by stad i the nearest AGN using i nter ferometry at infrared 

in polarised light in type 2 AGN (torus edge-on), most tests of wavelengths . The first objects to be successf ully analysed in this 



this "unified model" rely on the shape of the overall spectral en- way are the Seyfert j galaxy NGC 4 151 dSwainetal 



ergy distribution (SED) and statistics, i.e. the number of objects and the Seyfert 2 galaxy NGC 1068 (WittkowskTHi 



2003) 



. „ 2004; 

where the central source is obscured or visible. Until recently, a Jaflfe et al 2004; Ponce l et et al.ll2006h . In NGC 1068, which is 

direct spatially resolved assessment of the dust distribution that the brightest galaxy in the mid -infrared (MIR), indeed an ex- 

~~ ~ tended dust structure was resolved. In this paper we study the 

Based on observations collected at the European Southern . . • Ut » A/ -, XT • tU wtid *u r>- ■ 1 

, . „,„ . second brightest AGN in the MIR, the Circinus galaxy. 

Observatory, Chile, program numbers 073.A-9002(A), 060A-9224(A), & fe 3 

074.B-0213(A/B), 075.B-0215(A) and 077.B-0026(A). The Circinus galaxy is a highly inclined (~ 65°) SA(s)b 

Correspondence to: tristram@mpia.de galaxy harbouring a Seyfert type 2 active nucleus as well 
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as a n uclear starburst. At about 4 Mpc distance (1 arcsec ~ 
20 pc, Freema net al.l 1 1977b . the galaxy is among the nearest 
AGN and hence it is an ideal object to study the nuclear re- 
gion of active galaxies. The nucleus is heavily obscured by 
dust lanes in the plane of the galaxy, so that it is best ob- 
served in the infrared. The galaxy can be considered a pro- 
totype Sey fert 2 object due to the presence of narrow emis- 
sion lines (Oli va et al.ll 19941) . broad emission lines in polarised 
light d01ivaetal.lll998hT an ionisation cone traced by [Om], Ha 
and rSivnl dVeilleux & Bland-Hawthornlll997l: iMaiolino etall 
2000; IWilson et all 120001; iPrieto et alJ 120041) 7 an outflow ob- 
served in CO dCurranet all 1 199811 19991) . bipolar radio lobes 
( Elm outtie et al.ll 19981) and a Compton reflection component in 
X-ravs (lMattetal.lll996t ISmith & Wilson! 120011) . The ionisa- 
tion cone, the outflow and the radio lobes have a position an- 
gle of about -40°, which is roughly perpendicular to the disk 
of the galaxy, to the circum-nucle ar star forming rings a nd to 
the nuclear rings of molecular gas. iGreenhill et al.1 ([2003) have 
found an edge-on, warped "accretion" disk with a radius of 
'"outer ~ 0.40 pc which is traced by H2O masers. The almost 
Keplerian velocity curve yields an estimate of the mass of the 
central black hole to M BH < (1.7 ± 0.3) ■ 10 6 M . 

Modelling the SEP of the nucleus of t he Circinus galaxy was 
attempted bvlSiebenmorgen et all dl997l) . lRuiz et al.l d200l1) and 
ISchartmann et al.l (12005). The latter two came to the conclusion 
that the presumed toroidal dust distribution is relatively small, 
with an effective emission region in the MIR of less than 3 pc 
(i.e. < 0.15 arcsec). 

A detailed observational analysis in the mid-infrared was 
performed by Packham et al.l (120051) . The authors found ex- 
tended emission with 2 arcsec size at position angles of 8 1 
and 278°, which is consistent with the edges of the aforemen- 
tioned ionisation cones. However, their point spread function 
(PSF) analysis indicated that any extended nuclear emission in 
the mid-infrared has sizes less than 0.20 arcsec (4 pc) or is very 
weak. This showed again that single dish observations with cur- 
rent telescopes are not capable of resolving the nuclear dust dis- 
tribution even for the nearest AGN, such as Circinus. For this 
reason, interferometric observations of the Circinus nucleus in 
the mid-infrared with the Very Large Telescope Interferometer 
(VLTI) were performed. 

This paper is organised as follows. In Sect. |2]we present the 
details of the observations and of the data reduction. Sect. Ode- 
scribes the properties of the interferometric data and the results 
directly deducible from them. Sect. |4]treats the modelling of the 
data. The interpretation of our findings and implications thereof 
are discussed in Sect. [5] Our conclusions are given in Sect. [6] 



2. Observations and data reduction 

2.1. Instrument 

We obtained interferometric observations with the MID-infrared 
Interferometric instrument (MIDI) at the VLTI located on 
Cerro Paranal in northern Chile and operated by the European 
Southern Observatory (ESO). MIDI is a two beam Michelson 
type interferometer producing disperse d fringes in the N band 
in a wavelength range from 8 to 13 /jm dLeinert et all 2003). For 
our observations of Circinus, we combined the light of two of 
the 8.2 m unit telescopes (UTs) at any one time and used the low 
spectral resolution mode (A/6A « 30) by insertion of a NaCl 
prism into the light path after the beam combiner. 



2.2. Observations 

The data were acquired between February 2004 and May 2006 
during a total of six epochs. The first data set was obtained in 
February 2004 in science demonstration time (SDT). All follow- 
ing observations were carried out using guaranteed time obser- 
vations (GTO). Table \T\ gives a summary of the observations. 

The observing sequence repeated for the science target and 
calibrator star consists of the following steps: 

First, acquisition images of the two beams A and B coming 
from the unit telescopes are taken simultaneously. For this pur- 
pose, the standard chopping technique is used to suppress the 
background in the mid-infrared. For our observations, we used 
the default values for the chopping frequency (f - 2 Hz), chop- 
ping angle (a = 0) and chopping throw (5=15 arcsec) except 
for the observations on 2004 Jun 03, where a smaller throw of 
(5=10 arcsec was used. For the acquisition of Circinus, we used 
the SiC filter (11.8 ± 2.3 /mi) at longer wavelengths because of 
the rising spectrum of the source. For the acquisition of the cali- 
brator stars, which are in the Rayleigh-Jeans part of their spectra, 
we used the N8.7 filter (8.6+1 .5 /mi) at shorter wavelengths. The 
acquisition procedure was repeated until the overlap of the two 
beams was within one pixel accuracy, that is within ~ 86 mas. 

Then the beam combiner, the 0.6 arcsec wide slit and the 
NaCl prism are inserted into the light path leading to two spec- 
trally dispersed interferograms of opposite phase. No chopping 
is needed during the interferometric measurement as the back- 
ground is uncorrected and cancels out in the data reduction pro- 
cess. The fringe signal is searched by scanning the optical path 
difference (OPD) over a few mm at the expected position. This 
is achieved by adjusting the VLTI delay lines. At the same time, 
a piezo-driven mirror in MIDI performs a sawtooth scanning 
pattern with a scan length of 80/im for one of the delay lines. 
This modulation generates moving fringes from which the opti- 
cal path difference can be estimated. Once the zero optical path 
difference (0 OPD) is found, the so-called fringe track is started, 
during which the path difference is stabilised using the scanning 
of the piezo driven mirrors, while the VLTI delay lines com- 
pensate for siderial motion and atmospheric OPD shifts. There 
are two possibilities to track the fringes: the scans can pass over 
the white light fringe at OPD or they are performed next to 
it. We preferred to observe in "offset tracking" mode with the 
OPD stabilised at 50/iin from OPD, because this is of advan- 
tage for a coherent data analysis (see below). The observations 
on 2004 Jun 03 and 2005 Feb 21 and the first visibility point 
on 2005 Mar 01 were, however, observed tracking on the white 
light fringe. During these periods, the standard tracking software 
implemented at Paranal did not allow for using offset tracking. 

Finally, photometric data are recorded using only one tele- 
scope at a time in the otherwise identical optical set-up. This re- 
sults in two photometric data sets, one for each telescope. Again, 
chopping has to be used to suppress the background in the mid- 
infrared. 

Because our data were obtained over such a long pe- 
riod of time, the observation techniques evolved in the mean- 
time. Additionally, the instrument itself has improved consider- 
ably since the first observations were made in February 2004. 
This leads to a large variety of exposure times ranging from 
DIT = 12 ms to 20 ms. Similarly, the frame numbers range from 
NDIT = 5000 to 12000 for the fringe tracking and from NDIT = 
1500 to 4000 for the photometry (see Table [TJ. For the first 4 
visibility points, only a tip tilt correction for the i ncoming wave- 
fronts was possible using the STRAP units dBonaccini et al.l 
2000). A much better correction was obtained with the MACAO 



Table 1. Log of the observations, including the detector integration times (DIT), the number of frames (NDIT) for the fringe tracks and for the photometries as well as the lengths 
(BL) and the position angles (PA) of the projected baselines. Times, ambient values and baseline properties are for the start of the fringe track. 



Date and time 


Object 




DIT 


NDIT 


Airmass 


Seeing^ 


BL 


PA 


Associated calibrator and comments 


[UTC] 






[s] 


fringes 


phot£ 




[arcsec] 


[m] 


n 






2004 Feb 12: 


UT3 - UT2 






















06:55:11 


Circinus 


sciOl 


0.020 


12000 


(2000) 


1.473 


0.96 


43.49 


19.30 


calOl 


photometry unusable, used phot of sci02 


07:06:46 


Circinus 


sci02 


0.020 


12000 


2000 


1.449 


0.91 


43.36 


21.49 


cal02 


photometry used for sciOl and sci02 


07:36:34 


HD 120404 


calOl 


0.012 


12000 


2000 


1.456 


1.08 


41.39 


32.16 






08:27:38 


HD 120404 


cal02 


0.012 


12000 


2000 


1.416 


0.56 


40.33 


42.18 






2004 Jun 03: 


UT3 - UT2 






















05:50:45 


Circinus 


sci03 


0.012 


8000 


1500 


1.693 


0.68 


29.16 


92.99 


cal03 




06:42:24 


HD 120404 


cal03 


0.012 


8000 


1500 


2.095 


0.65 


25.40 


117.87 






07:58:04 


Circinus 


sci04 


0.012 


8000 


(1500) 


2.542 


0.73 


20.74 


129.68 


cal03 


phot B unusable, replaced by copy of phot A 


2005 Feb 21: 


UT2 - UT4 






















04:38:01 


HD 107446 


cal05 


0.012 


8000 


1500 


1.367 


0.83 


87.87 


51.88 






05:25:00 


Circinus 


sci05 


0.012 


8000 




1.617 


0.77 


87.38 


35.68 


cal05 


use photometry from 2005-03-01T04:38:30 


2005 Mar 01: 


UT3 - UT4 






















03:34:58 


HD 120404 


cal06 


0.012 


8000 


1500 


1.906 


0.95 


50.87 


39.76 




dsky = -2 (see Sect. 12. 3b 


04:06:04 


Circinus 


sci06 


0.012 


8000 


1500 


1.809 


0.75 


49.33 


44.58 


cal06 


dsky = -2 (see Sect. 12731 


04:38:30 


Circinus 


sci07 


0.012 


8000 


1500 


1.671 


0.72 


50.99 


53.61 


cal07 


photometry used for sci07 and sci08 


04:49:21 


Circinus 


sci08 


0.012 


8000 




1.631 


0.71 


51.57 


56.54 


cal07 


no photometry observed, used phot of sci07 


05:12:31 


HD 120404 


cal07 


0.012 


8000 


1500 


1.581 


0.55 


54.82 


65.93 






06:58:52 


Circinus 


sci09 


0.012 


8000 


1500 


1.358 


0.74 


58.09 


88.34 


cal08 




07:28:53 


HD 120404 


cal08 


0.012 


8000 


1500 


1.412 


0.48 


59.77 


98.71 






09:21:08 


Circinus 


scilO 


0.012 


8000 


1500 


1.345 


0.64 


61.95 


120.00 


cal09 


photometry used for scilO and scil 1 


09:43:25 


HD 120404 


cal09 


0.012 


8000 


1500 


1.480 


0.78 


62.10 


129.42 






10:08:22 


Circinus 


scil 1 


0.012 


8000 




1.400 


1.25 


62.36 


130.50 


cal09 


no photometry observed, used phot of scilO 


2005 May 26: 


UT2 - UT3 






















22:57:46 


HD 120404 


callO 


0.012 


8000 


1500 


1.666 


0.60 


42.62 


11.21 






23:29:03 


Circinus 


scil 2 


0.012 


8000 


3000 


1.563 


0.63 


43.77 


12.79 


callO 


reset chopping in second half of phot B 


23:43:07 


Circinus 


scil 3 


0.012 


5000 


3000 


1.524 


0.54 


43.67 


15.47 


callO 




01:37:29 


Circinus 


scil4 


0.018 


5000 


3000 


1.340 


0.67 


41.98 


37.05 


calll 




02:02:11 


HD 120404 


calll 


0.018 


5000 


1500 


1.409 


1.52 


39.64 


47.39 






03:40:07 


HD 120404 


call 2 


0.018 


5000 


1500 


1.454 


0.85 


36.24 


67.11 






04:07:02 


Circinus 


scil 5 


0.018 


5000 


3000 


1.374 


0.75 


36.82 


65.28 


call2 




04:17:27 


Circinus 


scil 6 


0.018 


5000 


3000 


1.387 


0.98 


36.31 


67.31 


call 2 




04:28:22 


Circinus 


scil 7 


0.018 


5000 


3000 


1.404 


1.05 


35.77 


69.45 


call2 




2006 May 18: 


UT2 - UT3 






















06:16:21 


Circinus 


scil 8 


0.018 


8000 


4000 


1.560 


1.50 


31.69 


84.18 


call3 




06:42:24 


HD 120404 


call 3 


0.018 


8000 


4000 


1.788 


0.84 


28.96 


100.07 






07:09:48 


Circinus 


scil 9 


0.018 


8000 


4000 


1.754 


0.98 


28.21 


96.40 


call3 




07:36:18 


HD 120404 


call4 


0.018 


8000 


4000 


2.035 


0.78 


25.98 


114.66 




rejected due to very low correlated flux 


08:04:13 


Circinus 


sci20 


0.018 


8000 


4000 


2.049 


0.59 


24.51 


110.75 


call3 


photometry used for sci20 and sci21 


08:17:34 


Circinus 


sci21 


0.018 


8000 




2.142 


0.67 


23.61 


114.69 


call 3 


no photometry observed, used phot of sci20 



H 



c 



X 



" From the seeing monitor (DIMM). 
A missing value signifies that no photometry was observed. A value in brackets denotes that the photometry was unusable. 
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adaptive optics system (lArsenault et al.l 120031) since its imple- 
mentation in late 2004, i.e. starting with our observations in 
February 2005. For Circinus, both systems were operated with 
a separate guide star at 50 + 2 arcsec distance from the nucleus. 
For the observations of the calibrators, the wavefront sensing 
was performed using the stars themselves. The overall ambient 
conditions for all the observations were good to fair with seeing 
values from the site monitor varying from 0.5 to 1.8 arcsec. We 
had no problems with clouds or constraints due to wind, except 
for one observation. This observation ("sci05"), carried out in 
service mode on the UT2 - UT4 baseline, had to be aborted due 
to clouds building up over Paranal. Thus, no photometric data 
could be obtained for this measurement. 

2.3. Data reduction 

Not all measurements contain the full data set consisting of a 
fringe track and two photometric measurements, so that they 
could be reduced in a straightforward manner. This has three 
reasons: first of all, several measurements were performed at the 
end of the night or of the allocated time and only the fringe track 
could be obtained. Secondly, for some measurements the pho- 
tometries proved to be unusable or to be of very low quality. 
Finally, for the measurement on the UT2 - UT4 baseline, no 
photometry could be obtained at all due to the adverse weather 
conditions. To nevertheless make use of all the fringe tracks ob- 
tained, we completed the fragmentary data sets with photomet- 
ric data from the other measurements preferably using measure- 
ments as close as possible in time. The details are included in the 
last column of Table [1] 

The data reduction was entirely performed by coherent visi- 
bility estimation using the software package EWS (Expert Work 
Station), written by Walter JafFe. The alternative software pack- 
age MIDI Interactive Analysis (MIA), which is based on a power 
spectrum analysis, was only used to cross-check the results for a 
selected number of observations: in these cases the two methods 
gave consistent results^. A detailed description of the coherent 
method can be found in |Jaffe|([2004). Here, we will concentrate 
only on the essential steps and highlight our specific configura- 
tion during the data reduction: 

First, a fixed mask is applied to the two opposite phased 
interferograms, which are then collapsed to form two one- 
dimensional interferograms. These are then subtracted with the 
result of doubling the signal while removing most of the back- 
ground. As the fringes vary rapidly in time due to the scanning 
of the MIDI internal mirrors, a high-pass filter in temporal di- 
rection with a width of, in our case, 10 frames (smooth =10) 
is applied to further remove the background. Additionally, the 
average of the entire spectrum is subtracted (dave = 1), assum- 
ing a modulation of the fringe signal in frequency space. This 
is indeed the case, when the OPD is not zero, hence our pref- 
erence to track next to the white light fringe. In the following 
two steps, the known instrumental OPD and the atmospheric de- 
lay are removed. The latter is estimated from the group delay 
(g smooth = 10). Additionally a constant phase shift induced by 
the varying index of refraction of water vapour is estimated and 
also removed. After rejecting frames with too low or too high 
optical path difference, the remaining data is averaged coher- 
ently yielding the raw correlated flux, rawcor, and the differen- 
tial phase, that is the relative variations of the phase with wave- 
length. 



1 An online manual of both MIA and EWS is available at 
http : / /www . strw . leidenuniv . nl/~koehler/MIA+EWS-Manual/ 



To determine the raw total flux of the target, first standard 
data reduction methods for chopped data are applied to the pho- 
tometric data: the sky frames are subtracted from the object 
frames. As this does not lead to a satisfactory sky subtraction, 
two stripes on the detector frame, one above the target position 
and one below, are used to subtract an additional wavelength de- 
pendent sky estimate. Only for one data set an adjustment of 
the position of these stripes using the option dsky = -2 was 
necessary (see the respective comment in Table [TJ. The data 
are then compressed and averaged to a one-dimensional pho- 
tometry, one for every telescope and side of the beam split- 
ter ("channel"). Two sets of total fluxes are extracted: One is 
needed for the determination of the visibility, and hence the 
same mask as for the interferometric data is applied. To ac- 
count for the possible inequality of the two beams from the two 
telescopes, this total flux is calculated as the geometric mean: 
rawtotg = V^i • B\ + y]A 2 ■ B 2 , with A\ the photo metry of 
beam A (from the first telescope) in channel 1, B\ the photom- 
etry of beam B (from the second telescope) in channel 1, as 
well as A 2 and B 2 the corresponding beams in channel 2. The 
second total flux is determined without applying a mask and 
by taking the arithmetic mean of the individual measurements: 
rawtota = (Aj + A 2 + B\ + B 2 )/4. The raw visibility is defined 
as rawvis = rawcor /rawtotg. The same data reduction is applied 
to both the science target and to the calibrator star leading to 
the raw correlated flux scirawcor and calrawcor respectively, as 
well as the raw total fluxes scirawtotg, calrawtotg, scirawtota, 
calrawtota and the raw visibilities scirawvis and calrawvis. 

Finally, the raw fluxes and visibilities of the science target 
are calibrated using the data of the calibrator star, to account 
for instrumental visibility losses and to determine the absolute 
flux values. For the calibration, we exclusively used the calibra- 
tor HD 120404 associated with Circinus, with the exception of 
one case, the long baseline UT2 - UT4, where HD 120404 was 
not observed at all. In this case, we used HD 107446 instead. 
The quality of the calibrators were checked by comparison to 
all the other calibrators observed during the same night. We de- 
cided to reject only one calibrator measurement because of its 
low correlated flux compared to the other measurements during 
the night: "call4" observed on 2006 May 18 at 07:36: 18. We cor- 
rected the Rayleigh- Jeans approximation for the template spec- 
trum of the calibrator using a template spectrum, caltemplate, 
from a database provided by Roy van Boekel, (private comm.). 
These template spectra were obtained by fitting stellar templates 
to multiband photometry. 

To keep the influence of errors in the photometry as 
small as possible (especially considering the smaller num- 
ber of usable photometries with respect to the number of 
successful fringe tracks), we calculated the calibrated corre- 
lated flux of Circinus according to F cm (A) = scicalcor = 
scirawcor/calrawcor - caltemplate using the option nophot = 1 
in EWS. The calibrated total flux is given by F tot (A) = 
scicaltota = scirawtota/calrawtota ■ caltemplate and the cali- 
brated visibility by V(A) = scicalvis = scirawvis /calrawvis = 
(scirawcor/scirawtotg) • (calrawtotg/calrawcor). Note that in 
this case V(A) * F COI (A)/F tot (*)- 

The entire data reduction and calibration was performed in 
the same manner for all visibility points without any further in- 
dividual adjustment of the parameters, in order not to introduce 
any biases. 
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Fig. 1. Coverage of the uv plane obtained for the Circinus galaxy with MIDI, the MID-infrared Interferometric instrument at the 
VLTI. The left panel shows the position of the observed fringe tracks in the uv plane for the two shorter baselines, UT2 - UT3 and 
UT3 - UT4. The points are coded according to the epoch of observation. The right panel shows all the 21 observed fringe tracks as 
well as the uv tracks of the baselines for the three combinations of unit telescopes that were used. 



2.4. uv coverage 

The interferometric data collected on Circinus constitute the 
most extensive interferometric data in the infrared or at shorter 
wavelengths of any extragalactic source to date. The coverage of 
the uv plane achieved by all of the MIDI observations is shown 
in Fig. Q] Each point represents the position of a fringe track 
leading to a visibility point. Because we measure the interfer- 
ometric signal of a real function (namely the brightness distri- 
bution) which has a symmetric Fourier transform, there are two 
points for each measurement, symmetrical to the centre. The dif- 
ferent colours denote the different epochs of observation. We ob- 
served 21 fringe tracks and 31 usable photometries (note this is 
not the double of the fringe track number as it should optimally 
be; there are only 15 usable A photometries and 16 usable B pho- 
tometries) for Circinus as well as 12 fringe tracks plus 24 pho- 
tometries for HD 120404, the calibrator star. From these data, we 
were able to reconstruct a total of 21 visibility points. 

The distribution of the measurements in the uv plane shown 
in Fig. Q] has concentrations at several locations, that is, two or 
three points at almost the same location. This was caused by con- 
secutive measurements of the interferometric signal. The differ- 
ences between such measurements can be considered as indica- 
tions for the accuracy of the measurements themselves. Several 
measurements also share the same photometry or calibrator data. 
In a strict sense, the true number of absolutely independent mea- 
surements is thus only on the order of 12, which is the number 
of independently measured calibrator data sets. 

The observations were performed with only three baseline 
configurations: UT2 - UT3, UT2 - UT4 and UT3 - UT4. This 
causes the visibility points to lie on arcs in the uv plane, which 
correspond to the classical uv tracks known from radio interfer- 
ometry. The uv tracks for these three baseline configurations are 
plotted as continuous lines in the right panel of Fig. Q] The uv 



coverage is far from optimal, as there are several regions lacking 
measurements: at a position angle of -20°, no measurements 
were obtained at all and at a position angle of +20° only one 
baseline length (BL = 40 m) was observed. The void regions 
therefore need to be filled in by future observations. In fact, the 
current coverage is unsuitable to attempt any image reconstruc- 
tion as commonly done with radio interferometry data. There, a 
flux distribution can be directly determined from the data by an 
inverse Fourier transform and a prescription, such as the "clean 
algorithm". We tried applying such an algorithm to our data. 
However, the form of the PSF (the "beam") corresponding to 
our uv coverage is very bad and shows strong periodic patterns, 
thereby dominating the outcome of the reconstruction process. 
The results of such an attempt are poor and hard to interpret, 
reflecting only general properties of the source. These general 
properties are much better constrained by models (see Sect. [4]). 

The differential phase (see Sect. 12.31 ) is less than 20° for al- 
most all baselines. At the shortest baselines, we find some ev- 
idence for a phase shift and there is also weak evidence for 
phase shifts at positions where the visibilities show unexpected 
behaviour (especially for sudden downturns at the edge of the 
spectal range covered). Nevertheless, we find no evidence for a 
major asymmetry in the source. We therefore postpone a detailed 
analysis of the differential phase information and concentrate in 
this paper on the interpretation of the visibility amplitude only. 

3. Results 

3.1. Total flux 

Figure [2] shows the spectrum of the calibrated total flux, F tot (A), 
of the Circinus nucleus. In several of the individual measure- 
ments (light grey lines), enhanced noisiness is apparent between 
9.5 and 10.0 yum, which is due to ozone absorption in the Earth's 
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Fig. 2. Spectrum of the total flux of the Circinus nucleus in the N band. The light grey curves are the 15 individual measurements 
of the total flux; error bars were omitted for the sake of clarity. The thick blue curve is the arith metic mean o f all th e individual 
measurements, the error bars indicating the reduced standard deviation. The spectrum measured bv lRoche et alj d!991l) is shown in 
green for comparison. The yellow curve shows the fit of the physical model to the total flux only (fit 1), while the red curve shows 
the fit to the full data set consisting of the total flux and the correlated fluxes (fit 2). The dashed red curve is the emission for the 
latter fit, without including the extinction by dust. For details on the modelling, see Sect. 14.21 



atmosphere. This is aggravated by the fact that the ozone feature 
coincides with the minimum of the flux in the intrinsic spec- 
trum of Circinus. The individual measurements are all consis- 
tent with the average (thick, blue line), when the associated in- 
dividual errors (not shown in Fig.|2| are considered. The curves 
with the largest deviations from the average come from mea- 
surements in 2004, which also have larger individual errors. The 
variations are mainly caused by imperfect background subtrac- 
tion for early data or observations at high airmasses. We clearly 
see an increase in the data quality after the first measurements. 
The edges of the N band are located at 8.2 /mi and 13.0//m. 
Towards shorter and longer wavelengths, the atmosphere has a 
high opacity due to water absorption. Given the large number 
of measurements, the averaged total flux is nevertheless usable 
beyond these boundaries. 

The broad absorption feature dominating most of the spec- 
trum is due to silicate absorption. We see no evidence for any 
line emission or emission of Polycyclic Aromatic Hydrocarbons 
(PAH ), as observ ed at larger distances from the nucleus 
dRoche et alj |2006) or in l arger apertures. The nuclear spectrum 
measured by Roche et ail agrees with our measurement of the 
total flux within 20 % (see Fig. [2j». The differences are probably 
due to the higher spectral re solution (R ~ 125) and the different 
orientation of the s lit in the Roche et alj d ata. The larger aper- 
ture data from ISO dMoorwood et alj|19 96) and Spitzer (unpub- 
lished) lie ~ 40 % above the MIDI total flux as they include the 
flux from the circumnuclear starburst. 



3.2. Interferometric data 

For all baselines observed, the visibility is well below the point- 
source value of V — 1, indicating that the emission region is well 



resolved with our interferometric resolution of A/2B < 40 mas. 
As an example, we show in Fig. [3] the correlated fluxes F C0I (A) 
and visibilities V(A) obtained on 2005 Mar 01 . During this night, 
observations were performed using the baseline UT3 - UT4. The 
baseline length only increased from 50 to 60 m, while the posi- 
tion angle underwent a major change from 44° to 1 3 1 °. As in the 
total flux, the broad dip at 10/mi is due to silicate absorption. 
The quality of the data at the edges of the wavelength range and 
in a window between 9.5 and 10.0/mi (hatched region) deterio- 
rates significantly due to atmospheric water and ozone absorp- 
tion. We therefore only trust the spectra within these boundaries, 
that is, from 8.2 /mi to 9.5 /mi and from 10.0 /mi to 13.0 /mi. 

A clear increase of the correlated flux and of the visibility 
can be observed with increasing position angle. This increase is 
most pronounced in the correlated fluxes longward of 1 1 /mi. We 
directly interpret this as an angle-dependent change of size of the 
emitting source: at position angles of 44° to 57° the correlated 
flux and the visibility are low, i.e. the source is more resolved and 
hence the emission extended, while at angles of 120° and 131° 
the correlated flux and visibility are higher, i.e. the source is less 
resolved and consequently the emission must be less extended. 
The angle-dependent size changes must be even stronger than 
Fig. [3] insinuates, as there is a small increase in baseline length 
towards larger position angles, that would lead to a decrease of 
the correlated flux and visibility for a round and smooth source. 
This is a direct and completely model independent evidence for 
an elongated dust structure. This elongation is perpendicular to 
the ionisation cone and the outflow at PA = -40° in the Circinus 
galaxy. 

To analyse the emissivity and temperature structure of the 
target either the visibility and the correlated flux can be used. 
Both approaches have advantages and disadvantages. The visi- 
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Fig. 3. Top panel: Spectra of the correlated fluxes, as observed during the night of 2005 Mar 01. Bottom panel: Visibilities for the 
same observations. During the course of these observations, the position angle (PA) changed by ~ 90°, while the projected baseline 
length (BL) stayed roughly the same. From the angular dependence of the correlated fluxes, direct implications on the morphology 
can be derived: the source observed is considerably more extended in the direction of PA ~ 50° than in the direction of PA ~ 130°. 
Areas affected by atmospheric absorption are hatched in grey. 



bility quantifies the degree of "resolvedness" of the source, that 
is, it contains geometrical information on the source. The corre- 
lated flux reflects only the emission from the unresolved part of 
the source, which per se does not contain any geometrical infor- 
mation. It only receives a geometrical meaning in comparison to 
the total flux or to correlated fluxes at different baselines. From 
our observations we know that the photometry is the most un- 
certain part of the measurement and we have fewer independent 
measurements of photometries than fringe tracks. As described 
in Sect. 12.31 the correlated fluxes are independent of errors in 
the photometry; the visibility, however, always contains the pho- 
tometry including its uncertainties. For the analysis and the mod- 
elling, we consider both the visibilities and the correlated fluxes 
as comparably useful and we will use both for our modelling of 
the data. Furthermore, we find that the essential inferences de- 
duced from either the visibilities or the correlated fluxes are the 
same; this is reassuring. 



4. Modelling 

To interpret the visibilities in terms of a flux distribution and 
to deduce physical properties of the source, the data need to 
be fitted by models, because the uv coverage is incomplete (see 
Sect. 12.4b and we lack absolute phase information. Our approach 
to model the data consists of several steps of increasing complex- 
ity, minimising any preconceptions about the source. The goal is 
to extract properties of the mid-infrared emission, such as size, 
elongation, depth of the silicate feature, etc., to sketch a picture 
of the nuclear dust distribution in the Circinus galaxy. 

The first step of modelling is aimed at estimating the char- 
acteristic size of the emission region using a simple, one dimen- 
sional model. From the discussion in Sect. 13 -2 1 it is clear that a 
one dimensional model is too simple to explain the MIR source 
in the nucleus of the Circinus galaxy correctly. More complex, 
two dimensional models are necessary. Our next step in mod- 
elling is hence a purely geometrically description of the source 
brightness distribution for each wavelength bin. This second, 
two dimensional model lead to the creation of a third, physi- 
cally motivated model which also attaches physical properties to 
the emission region. It turns out that the latter model compre- 
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hends the main results of the first two modelling steps. For this 
reason, only the physical model will be presented here while the 
descriptions of the one-dimensional and the purely geometrical 
model are found in App. [Aland App. IB"! 

All of the modelling was performed using the Interactive 
Data Language (IDL). 

4. 7 . Model selection 

For both the purely geometrical model as well as for the physi- 
cal model we explored a multitude of geometrical shapes in or- 
der to match the data. The shapes include point sources, unifom 
disks, elliptical disks, elliptical Gaussian distributions and com- 
binations thereof. We find that models including uniform disks 
produce results of similar y 1 values as Gaussian distributions. 
The brightness distribution of such models has sharp edges. The 
same is true for models including point sources {e.g. represent- 
ing the hot inner edges of the torus). Sharp edges have the great 
disadvantage that they produce oscillations in the uv plane: the 
uv plane is not uniform, but it has many minima and maxima in 
a quasi periodic pattern. Additionally, phase shifts occur at the 
minima in the uv plane. These are not observed. 

For this reason, we favour uniform models with no sharp 
edges to describe our data. This restricts our conclusions to 
more general statements about the properties of the mid-infrared 
source of Circinus. However, we are confident that these prop- 
erties are intrinsic to the source and neither an artefact of the 
specific model nor of the distribution of visibility points in the 
uv plane. To be able to make clearer statements about the de- 
tailed structure of the source, a more complete coverage of the 
uv plane is required and we are currently preparing to fill in the 
most obvious holes in the distribution of the visibility points, as 
far as this is possible with the baselines of the VLTI. 

4.2. 2D physical model 

The physical model consists of two concentric black body 
Gaussian emitters behind an extinction screen. For the absorp- 
tion we used the wa velength dependency of th e extinction of in- 
terstellar dust from Sch artmann et al.l (2005), but modified the 
profile of the silicate absorption featur e in the wave l ength range 
from 8.0 to 12.7 yum using data from Kem per et alj (|2004). The 
latter determined the feature profile towards the galactic centre 
and we thus expect their profile to be our best guess for the line 
of sight towards the nucleus of the Circinus galaxy. To obtain our 
template absorption profile t(A), the entire optical depth curve 
was normalised to 1 at the maximum absorption depth of the 
silicate profile: r(9.7yum) = 1. 

Our two component model has the following functional 
form: 



F(A,a,5) = fi -Gx{A,a,5)-&- TvTW 
+ f2-G 2 (A,a,S)-e-^ TW 



(1) 



where T\, t 2 are the individual optical depths of the absorption 
screens. The scaling factors f\ and f 2 are factors to scale the 
intensity of the black body Gaussian distributions G\ and G 2 . 
The Gaussian distributions in turn are defined by 



G n (A, a, 6) = exp 



- 4 In 2 • ^ u-cos <A+^sin0 'j 2 

- 4 In 2 ■ ^±El£dLE2l£^ 2 



(2) 



Table 2. Best fit parameters for the two component black body 
Gaussian models. The first column is for a fit of the total flux 
only (fit 1), the second column for a fit to both the total flux 
and the correlated fluxes (fit 2). The values in brackets were not 
fitted, but held fixed. 



parameter 


fit 1 


fit 2 


Ai [mas] 


23.6 


21.1 


T\ 


(0.25) 


0.21 


Tl 


1.17 


1.18 


r, [K] 


479.6 


333.7 


fx 


(1.00) 


1.00 


A2 [mas] 


279.0 


96.7 


ri 


(1.00) 


0.97 


T 2 


2.85 


2.22 


T 2 [K] 


144.4 


298.4 


h 


(1.00) 


0.20 




(60.0) 


60.9 


X 1 IN free 


0.56 


36.86 



for each of the two components n = 1,2. The parameters of the 
Gaussians (r\, r 2 for the axis ratios, Ai and A2 for the FWHM 
and <p for the orientation) are the same as for the purely geomet- 
rical model (see App. |B), except that they are not wavelength 
dependent. is the black body intensity which depends on the 
temperature T„. This model has 11 free parameters: the FWHM 
A„, the axis ratio r„, the temperature T„, the flux normalisation /„ 
and the optical depth r„ for each of the two components n = 1,2. 
The common position angle <p is the eleventh parametefl Using 
these parameters, the flux distributions at several wavelengths 
were generated and Fourier transformed. From the resulting 
distribution in the uv plane, the visibilities corresponding to 
the observed baselines, Vf° & {A) = \Vx(u u v;)l, were extracted. 
Additionally, the total flux was determined by integrating over 
the entire emission region: F™° d (A) = J J F(A, a, 6) da d<S. The 
comparison to the observations was performed using both the 
correlated fluxes and the total flux, because here the precise val- 
ues should match and not only the flux ratios. The correlated 
fluxes for the model were obtained by multiplying the mod- 
elled visibilities with the total flux of the model, i.e. F mod (/l) = 
V™ od (^) ■ F™ d (A). These could then be compared to the data. 
The model parameters were optimised using the Levenberg- 
Marquardt least-squares minimisation. The % 2 to be minimised 
was: 

21 ( robs /" i \ i7mod / i \ \ 2 



X 



=zz 



( F° hs (A) - F™ oi (A) 

(T p-obs (A) 



(3) 



where the sum over i is the summation over the 22 observed 
baseline orientations (21 baselines plus the total flux as a base- 
line with a baseline length of BL = m). The cr F? b S (/l) are the 
errors of the correlated fluxes. 

Two fits were performed: first the model was only optimised 
to fit the total flux (in the following "fit 1"). In this case, several 
parameters (r\, r 2 , f\, f 2 and 0) were held fixed as they were de- 
generate. Then the model was fitted to the total data set including 
both the total flux and the 21 correlated fluxes ("fit 2"). In case 
of the fit to the total flux only, there is hence no summation over 
i in Eq.|3] as we simply had to use i = 0. The result of the fitting 
process is summarised in Table [2] 



Fbb(T„, A) 



1 Decoupling tp for the two components leads the fits to converge on 
a few dominant groups of visibility points only. This is the reason for 
keeping a single orientation angle for both components. 



Tristram et al.: The nuclear dust torus in the Circinus galaxy 



9 




Fig. 4. Comparison of the correlated flux predicted by the physical model (fit 2) with the observations for two distinct wavelengths: 
(a) 8.5 and (b) 12.5 yum. The measured fluxes are plotted with asterisks; the continuous lines are the modelled fluxes. Note the 
different flux ranges for the different wavelengths. 



A total of 21 wavelength bins were used for the fitting, so 
that fit 1 has A^f ree ,i =21-1 — 11 = 10 degrees of freedom, while 
fit 2 has -/Vfree.2 = 21-22 -11 = 451. The resulting total fluxes 
for the two models are plotted in Fig. [2] in orange (fit 1) and red 
(fit 2) on top of the measured total flux. The observed correlated 
fluxes are compared to fit 2 in Fig. [4] at the two wavelengths 
already previously featured: 8.5 and 12.5 fim. A comparison of 
fit 2 with all of the data is given in App.lCl 

The two components of this model have similar properties 
to those of the purely geometrical model (App. |B): there is a 
smaller elongated component and a larger, nearly round com- 
ponent. A sketch of the result of fit 2 is depicted in Fig. [5] 
and the flux distribution at 1 1 fim can be seen in Fig. [6] The 
smaller component has a FWHM of Ai =21 mas, which corre- 
sponds to 0.4 pc at the distance of Circinus. It is highly flattened 
(n = 0.21) and has a temperature of T\ = 334 K. The distri- 
bution is found to be an optically thick black body ( f\ = 1). 
The second component is significantly larger, as a large part of 
the observed flux (~ 90%) is resolved with our interferomet- 
ric set-up. It has a FWHM of A2 = 97 mas, which corresponds 
to 2.0 pc, only a very small ellipticity and a temperature of less 
than 300 K. The emissivity of the black body radiation from this 
component is scaled by a factor of = 0.20. This low scaling 
factor is driven by the rather extended flux distribution, which 
is necessary to explain the low visibilities, and the requirement 
not to overpredict the measured total flux. Note that the scal- 
ing factor is not to be interpreted as a geometrical covering fac- 
tor. Together with the charactersitic temperature T2, it rather 
is a formal description of a significantly more complex multi- 
temperature system. This is underlined by the fact that there is a 
degeneracy between the temperature and the scaling factor in the 
fit. When fitting only the total flux, we see a larger temperature 
difference between the two components than for the full fit. In 
both cases, fit 1 and 2, the temperatures are the least constrained 
parameters. We therefore consider the 330 K for the inner com- 
ponent as a lower limit for the highest dust temperature and the 
300 K as an upper limit for the cool component. Our data rules 
out any significant contribution from a truly hot component with 
temperatures T > 1000K, i.e. close to the sublimat i on tem pera- 
ture of the dust. This is consistent with lPrieto et alj d2004). who 
also found no evidence for hot dust. The conclusion drawn from 
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Fig. 5. Sketch of the physical model (fit 2): a highly elongated 
warm emission region with a temperature of ~ 330 K (yellow) 
is surrounded by an extended, almost round and slightly cooler 
emission region with a temperature of ~ 300 K (brown). In the 
centre of t he sketch, the location of the H2O maser emitters in a 
disk from iGreenhill et ail (120031) is overplotted: the blue line to 
the north-east traces the receeding masers and the red line to the 
south-west the approaching masers. The dashed light blue line 
traces the edge of the observed ionisation cone, the dotted line is 
the cone axis. 



the growth of the size of the emitter with wavelength (Sect. lAl . 
namely that the temperature of the emitter decreases with in- 
creasing distance to the nucleus, is confirmed by the lower tem- 
perature of our larger component. The small component only 
contributes a minor fraction to the total flux. This can be seen 
in the top left graph of Fig. lCTI the flux contribution by the ex- 
tended component is delineated by a dashed line and that by the 
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Fig. 6. Flux distribution of the physical model (fit 2) at 1 1 /mi in Fig. 7. Flux distribution of the physical model with added 



logarithmic scaling. 



dumpiness for the extended component at 1 1 /mi in logarithmic 
scaling. 



compact disk component by a dotted line . In contrast to the ob - 
servations of NGC 1068 dJaffe et al.ll2004t iPoncelet et alj|2006h . 
the correlated fluxes in the Circinus galaxy do not show a rela- 
tive increase of the flux to shorter wavelengths with respect to 
the flux at longer wavelengths. This shows that there is only a 
moderate increase of dust temperature towards the centre. 

The model used here is a parametrised global model and 
does not match all aspects of the data. The reduced x 1 of 37 
indicates that deviations well above the 5<x-level are commonly 
found. The aim of our model was to extract the overall properties 
of the emitting source. The general behaviour of the correlated 
fluxes for the different baseline orientations can be well repro- 
duced, but not the details. Looking at Fig. 2J the model fluxes 
only show a single dip in the range of position angles traced, due 
to the symmetry of the object when rotated by 180°. The mea- 
surements exhibit considerably more rapid angular variations, 
rejecting a smooth axisymmetric model. The dispersed fluxes 
(Fig. IC.ll ) disagree most strongly at the edges of the spectrum, 
at 8 /mi, where a downturn in several of our correlated fluxes is 
seen, or at 13 //m, where the model deviates in both directions 
with respect to the data. As there are also small changes in the 
differential phase in these regions, we expect a further substruc- 
ture of the source to be responsible for this unusual behaviour. 

To verify this possibility, we undertook an additional mod- 
elling experiment: we tried to reproduce the low scaling factor 
of the larger component through dumpiness: instead of multi- 
plying F\ with the factor fa, an uneven surface C(a, 6) was used. 
C(a, 6) can be interpreted as a screen simulating dumpiness. To 
make the calculation consistent, the average value of C has to 
reproduce the scaling factor: (C(a, 6)) = fj . The resulting flux 
distribution for such a clumpy model is shown in Fig. [7] In such 
a simple experiment, it is possible to find clump distributions 
which trace the correlated flux much better (see for example 
Fig- The reduced^ 2 drops to values of x 2 IN f ree < 25. This 
supports the assumption that such deviations are generated by a 
substructure of the brightness distribution. 



5. Discussion 

The modelling of our interferometric observations reveals the 
presence of two components of dust emission at the very cen- 
tre of the Circinus galaxy: an extended (r\ = 1.0 pc), warm 
(J2 ~ 300 K) and fairly round emission region showing strong 
silicate absorption and a much smaller (ri = 0.2 pc), but only 
slightly warmer (T\ ~ 330 K) component that is highly flattened 
and shows a less pronounced silicate absorption feature. We in- 
terpret these two components as signs for a geometrically thick 
"torus" of warm dust surrounding a warmer, disk-like structure 
(see Fig. [9]). The smaller and denser disk is seen at a high inclina- 
tion and it partly exhibits the silicate band in emission. It is sur- 
rounded by the larger, less dense and most likely clumpy torus 
component which gives rise to strong silicate absorption. The 
silicate absorption profiles towards both the disk and the torus 
component are consistent with the interstellar absorption profile 
observed towards the nucleus of our own galaxy. This suggests 
that the dust composition is the same. The result is interesting 
compared to that found for NGC 1068, where the dust absorp- 
tion towards the nucleus clearly does not agree w i th the standard 
intergalactic absorption profile dJaffe et alj|2004t IPoncelet et al.l 
I2006h . 

To understand the structure of the nuclear dust distribution 
in AGN, we developed hy drodynamical mod els of dusty tori (for 
a detailed description see Schartmann 2007). These models in- 
deed also show a disk-like structure in the inner region which 
is surrounded by a geometrically thick "torus". The torus is not 
continuous but has a filamentary structure (rather than clumpy) 
and the disk is turbulent. 

In the following, we will develop a picture of the central 
parsec region in Circinus which aims to explain the results of 
the MIDI measurements in the context of the hydrodynamical 
model. In addition, the picture needs to be consistent with the 
previously known properties of this galactic nucleus. 
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Fig. 8. Correlated flux of the physical model with simulated dumpiness for the extended component: (a) 8.5 and (b) 12.5 fim. The 
measured fluxes are plotted with asterisks; the continuous lines are the modelled fluxes. The dumpiness introduces significant 
changes in the correlated fluxes compared to the smooth model (see Fig. |4j». From among 1000 clumpy models we found at least 
one (shown here), which matches the data well. This shows that dumpiness of a certain geometry can greatly improve the fit to the 
data. 




Fig. 9. Cut through the dusty torus in the central parsec-sized 
region of the Circinus galaxy derived from the interferometric 
observations in the MIR. See text for details. 

5.1. Extinction 

To interpret the nuclear emission from Circinus and the silicate 
absorption in the torus, we need to estimate the foreground ex- 
tinction. Radiation from the nuclear region of Circinus suffers 
extinction from two main foreground absorbers: our own galaxy 
(Circinus is located at a galactic latitude of b — -4°) and the 
intrinsic absorption in the galactic disk of the Circinus galaxy. 

By obser ving the reddening o f several stars in the vicinity 
of Circinus. iFreeman et al.l d 1977I) estimate the visual extinction 
due to our own galaxy to be only Ay - 1.50 ± 0.15 mag, as 
Circinus lies in a window of low galactic extinction. The extinc- 
tion towards the nuclear region within Circinus itself is much 



less constrained and estimates from near-infrared colors vary 
from Ay — 6 mag for a foreground screen to Ay > 20 mag for a 
mix o f dust and stars (iPrieto et al.ll2004t iMueller Sanchez et al.l 
120061) . In the latter case, the colours saturate so that higher ex- 
tinctions cannot be ruled out. The near-infrared observations 
which underlie this estimate probe a region on scales of less than 
10 pc dist ance from the nucleus. Given the presence of a nuclear 
starburst (Muel ler Sanchez et al. 1 12009) . a mix of dust and stars 
appears to be more realistic and we adopt a foreground extinc- 
tion, between us and ~ 2 pc from the nucleus, of Ay ~ 20 mag. 

From the optical depth of the silicate feature seen in the ex- 
tended torus component, we can estimate the extinction towards 
those regions of the torus where the warm dust emission arises. 
Our value of T2 ~ 2.2 for the optical dep th is at the lower end 
of the values found bv lRoche et al.l (120061) . who cite a range of 
2.2 < T S i] < 3.5, depending on where they extract the spec- 
trum along their slit. Given the low spectral resolution of MIDI 
the minimum of the absorption trough may be washed out, lead- 
ing to lower optical depths. We therefore consider the two val- 
ues to be in good agreement, although our observations rule out 
T s ii > 3.0 for the nucleus. In the following, we will assume 
T s ii = 2.5 towards the larger (r = 1 pc) torus component. 

We use the opacity scaling of Schartmann et al.l (120051) . to 
convert the opacity measured in the 9.5 /mi silicate feature to 
an equivalent visual extinction, Ty = 12.2 ■ r s ji. Hence, the ab- 
sorption in the visual is Ay = 1.09 ■ Ty mag = 13.3 ■ T s n mag. 
For T s ii = 2.5, we obtain Ay = 33 mag, which means that there 
are an additional 1 1 mag of extinction towards the torus region 
compared to the limit derived from the near-infrared colours. 
Most likely, the extinction towards the accretion disk itself is 
even higher. 

In contrast to the MIDI obse rvations of NGC 1068 
dJaffe et alJl200H iPoncelet et all2006h . we do not observe a fur- 
ther increase of the silicate absorption feature when zooming 
into the nucleus with our interferometric observations. Instead, 
the absorption feature appears to be less pronounced in the corre- 
lated fluxes, which can also be inferred directly from the convex 
shape of the visibilities (e.g. from Fig. |3J. In general, the sili- 
cate feature is expected to appear in emission for lines of sight, 
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where the (appropriately high) temperature decreases away from 
the observer. This is the case for the directly visible inner rim of 
the torus and, therefore, expected and observed for Seyfert 1 ob- 
jects. In the opposite case (cold dust along the line of sight), the 
silicate feature will appear in absorption, as seen in Seyfert 2 
objects. A decrease of the silicate feature in absorption therefore 
means that more and more of the Seyfert 1 type silicate emis- 
sion is able to shine through the torus body. These inferences are 
consistent with a hidden Seyfert 1 nucleus in Circinus, a fact also 
sugg ested by the observ ation of broad emission lines in polarised 
light dOliva et al.|[l998h . 

In summary, we conclude that the obscuration on our line 
of sight towards the inner parsec of Circinus is composed 
of the following three extinction components: extinction by 
our Milky Way, Ay(MW) = 1.5 mag, by the foreground dust 
within Circinus, Ay(disk) = 20 mag, and by the nuclear dust, 
Ay(nuclear) =11 mag. 

5.2. Energy budget 

The energy emitted by the source in the N band can be estimated 
by integrating the flux in the dereddened spectrum, giving a lu- 
minosity in the N band of L 7 .5-i3.5^m = 1.3-1O 9 L = 5.1-10 35 W. 

To estimate the total luminosity of the dust emission, we inte- 
grate under the flux of our physical model (fit 2) from v = to oo 
without applying any absorption, that is t\ — T2 = 0. Orientation 
effects only affect significantly the smaller disk component. A 
correction by assuming oblate spheroids seen edge-on leads to 
an increase of the flux by about 20 %. With this method, we ob- 
tain a total luminosity from the dust of Ljust = 4.1 • 10 9 L Q . This 
is four times more than the luminosity measured in the N band 
alone. Considering that we have not traced the contribution by 
dust at T < 150 K, we round our value up to Ld ust ~ 5 • 1O 9 L 
as our best estimate for the luminosity of the nuclear dust at 
r < 2 pc in the infrared. 

We can take this value to estimate the luminosity of the 
accretion disk, L acc , which presumably is responsible for heat- 
ing the dust. In doing so, we have to take into account that 
not all of the accretion luminosity is available for heating the 
torus. A significant fraction of L acc will escape along the di- 
rections of the ionisation cones. With an opening an gle of 90° 
dVeilleux & Bland-Hawthornll 19971: IWilson et al.ll2000h . the to- 
tal solid angle covered by the torus, as seen from the cen- 
tre, is 2 VZtt. This corresponds to about 0.7 of the full sphere. 
Assuming a radiation characteristic for an optically thick accre- 
tion disk which is proportional to cos#, about half of the ra- 
diation is emitted in direction of the dust in the torus. That is, 
only half of the energy from the accretion disk will be absorbed 
by the dust and re-radiated in the infrared. If we thus assume 
idust ~ 0.5 • Lacc, this leads to a luminosity of the central energy 
source of L acc ~ 1O 1O L . 

Consistent with our estimate, iMoorwood et al.l (1 1996b also 
report the luminosity of the nucleus to be ~ 1O 1O L . Their inte- 
gration under the UV bump (the unobscured radiation) yields the 
same luminosity a s we observed in the MIR: Luv = 5 ■ 1O 9 L . 
Oliv a et al.l ( 1 19991) cite a total ionising radiation of the nucleus on 
the order of Li on = 1 -4- lO lo L -sin i, where i is the inclination 
angle o f the ionisation cone w ith respect to our line of sight. For 
i > 60 dElmouttie et alJI 19981) . this translates into 1 - 7 • 10 10 L . 
Luminosities of the AGN in Circinus significantly higher than 
10 10 L can be ruled out, as this would imply a higher MIR flux 
than observed. The currently derived intrinsic 2 to lOkeV lu- 
minosities of the nuclear X-ray source are all below 5 • 1O 8 L , 



which is well below the energy needed to heat the dust. The dust 
is mainly heated by soft X-rays and UV radiation rather than by 
radiation in the 2 to lOkeV band alone. 

Finally, we can compare our luminosity estimate for the ac- 
cretion disk to the Eddington l uminosity. For a blac k hole with 
a mass of M BH = 1.7 ■ 10 6 M dGreenhill et al.ll2003l) . this lumi- 
nosity is L E = 5.6 ■ 1O 1O L . With our value of L acc ~ 1O 1O L , 
we find Circinus to accrete at 18 % of its Eddington luminosity. 

5.3. Temperature dependency and dumpiness 

Less than a parsec from this central energy source, the dust is 
most likely heated by the radiation from the central accretion 
disk alone and not by a homogeneous energy input, such as by a 
starburst. The discussion of the energy budget of the nucleus has 
shown that no further energy source is necessary. A heating by 
a nuclear energy source is also suggested by the decrease of the 
dust temperature to the outside. 

The temperature decrease we observe is, however, very shal- 
low and the temperature of our large component is relatively 
high considering the distance of its dust to the heating source: 
the dust has a temperature of ~ 300 K out to distances of 1 pc 
from the nucleus. To investigate the implications of this finding 
on the properties of the dust distribution, we compare the dust 
temperature in Circinus to the temperature expected for a direct 
illumination of the dust and to the temperature inferred from our 
radiative transfer modelling of continuous and clumpy AGN tori. 
The various temperature dependencies are shown in Fig. [10] and 
we will discuss them one by one in the following. 

Assuming a naive radiation equilibrium, the temperature of 
a dust particle at the distance r from the illuminating source with 
the luminosity L acc is given by 

Taking into account the absorption and emission characteristics 
of the dust grains leads to a more realistic description of the 
temperature dependency (lBarvainislll987l) . For emission in the 
infrared, the temperature dependency then reads: 

The two temperature dependencies T^{r) and Tg(r), calculated 
for L acc = 10 10 L , are plotted Fig.[10]by a green dashed line and 
a green continuous line respectively (in the printed version the 
green lines appear grey). The realistic case leads to significantly 
higher temperatures than the naive case, as the individual dust 
particles reradiate the energy less effectively than a black body. 
Here, we have not taken into account the anisotropic radiation 
characteristic of the central energy source, the accretion disk, 
which emits less energy in direction of the dust distribution than 
in direction of the ionisation cones (see Sect. 15.21 ). This leads to 
a lower temperature of the dust near the plane of the disk, so 
that Ts(r) should be considered as an upper limit for the dust 
temperature at a certain radius. 

In the case of a continuous dust distribution, the tempera- 
ture is expected to drop steeply outside the inner funnel of the 
torus, where the dust is directly illuminated by the central source. 
The dust at larger radii is shielded from the nuclear heating and 
by consequence at lower temperatures. For a clumpy torus, the 
temperature dependency should be somewhat intermediate be- 
tween the continuous case and the radiation equilibrium for to- 
tally unobscured lines of sight, that is, T%(r). The behaviour is 
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Fig. 10. Radial temperature gradients in the central dust distribution. The radial temperature distribution of our simple "physical" 
model (thick black continuous lines) is compared to the temperature distributions for a clumpy (hatched area) and continuous torus 
(light blue / grey filled area) as well as to the temperature dependency of dust directly illuminated by the central energy source 
(green / grey lines). The point where the sublimation temperature of the dust, Tb(t) - 1500 K, is reached is marked by a green / 
grey circle. The lower part of the graph shows the intensity of the Gaussian components (depending on the scaling factor /). The 
HWHM of the Gaussian components are marked by circles. The scales probed by our interferometric set-up lie within the dotted 
vertical lines. 



reproduced by our radiative transfer modelling of AGN tori: the 
shaded areas show the temperature ranges (from mean to max- 
imum temperature) of the individual cells in a continuous and 
a clumpy torus model (for a detailed de scription of the mod- 
elling we refer to lSchartmann et aT1l2007l) . For this purpose, the 
radii were rescaled so that the inner radii of the torus, where 
the sublimation temperature is reached, coincide with the subli- 
mation radius predicted by Tgir). As expected, the continuous 
model (filled, light blue / grey area) has a very strong drop in the 
temperature behind the inner edge of the dust distribution, while 
the clumpy model (dark blue, hatched area) has a much wider 
temperature range, including lower temperature clouds at small 
distances from the nucleus and higher temperatures at larger dis- 
tances. 

Finally, the thick black continuous lines trace the temper- 
atures of the two dust components in Circinus. The lines are 
drawn out to the half width at half maximum (HWHM) of the 
respective component. From the figure it is obvious that the tem- 
peratures seen by MIDI disagree with a continuous dust distri- 
bution as the temperatures lie significantly above the distribu- 
tion and no sharp temperature rise to the centre is seen. The 
high temperatures at large distances from the nucleus can only 
be achieved when clouds also at larger distances have free lines 



of sight to the heating source. As the temperatures lie within the 
temperature range for a clumpy dust distribution, we conclude 
that the dust most likely is distributed in a clumpy medium. All 
our temperatures lie below those of Ts(r), as expected for partial 
shadowing of the clouds. 

There are several further indications for a "clumpy" or 
patchy dust distribution in our data: 

The low emissivity filling factor of the torus-like component 
in our physical model is required by the data and best explained 
by an inhomogeneity of the emitting dust. Patchy emission ap- 
pears the only way for thermal emission of dust at T ~ 300 K 
to be as extended as observed here, while at the same time being 
consistent with the total flux of the source. 

The most compelling evidence, however, comes from the be- 
haviour of the correlated fluxes (a) with wavelength or (b) with 
changing baseline (for the sake of simplicity we will only refer 
to the correlated fluxes F cor in this paragraph, nevertheless the 
same statements are true for the behaviour of the visibilities V). 
The behaviour observed cannot be explained by a smooth bright- 
ness distribution. 

When considering the set of all visibility points observed at 
a single wavelength, i.e. F cor / \(PA), the need for dumpiness be- 
comes obvious when comparing Fig.|4]to Fig. [8] as was already 



14 



Tristram et al.: The nuclear dust torus in the Circinus galaxy 



discussed in Sect. 14.21 We have also compared the variations we 
observe in our data to those that can be see n in the correlated 
fluxes of our clumpy radiative transfer models (ISchartmann et al.l 
120071) : when we calculate the correlated fluxes of such a torus 
model for different position angles, F cor j(PA), we obtain varia- 
tions of up to 50 % of the flux, similar to those shown in Fig. [8] 
The variations become stronger for longer baselines, where sub- 
structures are better resolved. 

Similarly, the need for a clumpy structure becomes obvi- 
ous when considering the individual visibility points at different 
wavelengths, that is, the dependence of the correlated flux over 
the N band, F C or,.m(/l)- In this case, we observe wiggles in the 
correlated fluxes that can only be reproduced by a substructure, 
e.g. the downturn at A > 12 /mi for BL ~ 36 m, PA ~ 67° (see 
Fig. ED. 

An inhomogeneous dust distribution complicates the prob- 
lem of determining the morphology of the source significantly, 
when interferometric methods are used. With the current uv cov- 
erage, we are not able to say anything about the details of the 
clump distribution of the source. More data are needed to disen- 
tangle these and also to enable us to trace how the morphology 
changes with the wavelength. 

5.4. Orientation and geometry of the dust distribution 

Both direct inspection of the data (see Fig.O and our modelling 
reveal that the source is elongated in direction of PA ~ 60°. This 
effect is dominated by our observations with longer baselines 
(UT3 - UT4); the shorter baseline observations (UT2 - UT3) do 
not show a clear orientational preference. In our interpretation, 
this means that the inner component is more disk-like, while the 
outer component traces an almost spherical or very thick toroidal 
component. 

The projected orientation of our disk component is in very 
good agreement with the orientation of several components in 
this AGN which are directly related to the nuclear dust distri- 
bution: its axis coincides with the directi on of the ionisation 
cone, which point s towards PA = -44° (Maiolino e t al]|2"o"o"ot 
Wils on et al.l20 00) and the bidirectional out flow observed in CO 
pointing towards PA = 124° and PA = -56° dCurran et al.ll998l) . 
Our measurements are hence the first direct confirmation for the 
presence of a dust structure extended perpendicularly to the out- 
flow and the ionisation cone in an AGN without taking any pre- 
vious knowledge on the source into account. 

The nucleus of Circinus is special in another feature: it dis- 
plays a very strong emissio n of water vapour mas e rs. Th is emis- 
sion was first detected by Gardn er" fe Whiteoakl (Il982l) . A de- 
tailed study was performed bv lGreenhill et al](l2003l) : The VLBI 
observations of the 6 k,- 523 transition of H2O show an S shaped 
locus and a wider distribution of maser sources at the nucleus of 
the Circinus galaxy. The authors attribute the sources detected 
to two distinct populations: a warped, edge-on "accretion disk" 
and a wide angle outflow. The warped disk has an inner radius 
of (0.11 + 0.02) pc and extends out to ~ 0.4 pc. These are the 
only observations at the same or even better angular resolution 
as ours. An almost perfect alignment is found between the outer 
radius of the maser disk with a position angle of PA = 56 + 6° 
dGreenhill et al.ll2003l) and our disk component (c.f. Figs. [5]and 

Because of the high axis ratio of r\ = 0.21, our observations 
suggest that the disk component is nearly edge-on: i'disk > 60°. 
If it were a truly edge-on disk (z = 90°), it would then have to 
be warped or to have a thickness of h/r = r\ to appear as we 
observe it. When observed edge-on, the maser disk with a warp 



of 27° can easily produce the same extent as observed for our 
smaller component. Considering this and the good agreement of 
the position angles, we believe the maser emission originates in 
the dense dusty and molecular medium associated with our disk 
component 

From the present data, we cannot make any detailed state- 
ments concerning the shape of the larger torus component, as 
its proportions are not very well constrained. It must be signifi- 
cantly extended also in the direction of the ionisation cone and 
the measurements are consistent with an almost circular dust dis- 
tribution when seen in projection. When seen almost edge-on, 
this implies h/r ~ 1. The component can, however, not be en- 
tirely spherical but it must have an opening angle of 90° for the 
ionisation cones (see Fig.[9j>. From our side-on viewing position, 
we see a projection of this hollowed spherical dust structure, 
which leads to its almost circular appearance. More measure- 
ments at shorter baselines are necessary to provide better con- 
straints on the distribution of the more extended dust emission. 

5.5. Collimation of the ionisation cone and outflow 

In the context of the two components, namely the disk and the 
torus, it is interesting to address the question which component 
causes the collimation of the outflow and the ionisation cone. 
The latter has an opening angle o f 6\ = 90° and a sh arp bound- 
ary seen clearest in [Om] and YLa (Wil son et al.l20 00). The same 
opening angle is found for the e xtended X-ray emissio n, over- 
lapping with the ionisati on cone (Smith & Wilson 2001) as well 
as for the CO outflow dCurran et al.l 1 19991) . Finally, the high 
ionisation lines in d irection of the ionisation cone re ported by 
IMaiolino et ail (l2000h ([Sivi]) and iPrieto et all (120041) ([Si vn]) 
show an even smaller opening angle. Summarised, this means 
that at scales > 2 pc only 30 % of the solid angle are exposed 
directly to radiation from the nuclear engine. This cannot be 
explained by the cos 6 radiation characteristic of the accretion 
disk, especially consi dering the sharp bound aries of the cone. On 
scales less than 1 pc. iGreenhill et al] (12003) found a wide angle 
outflow of maser emitters. It is supposed to cover 80 to 90 % of 
the solid angle, corresponding to an opening angle of 9 m = 120°. 
This is consistent with a collimation of the outflow by a disk with 
a warp of 27°. These observations lead to the conclusion that the 
large torus component extending from ~ 0.2 to 1 .0 pc must be 
responsible for the further collimation of the wide angle outflow 
to the narrow one observed at larger scales (see Fig.|9j. 

Further support for this picture comes from the observa- 
tion of molecular hy drogen gas just below the ion isation cone 
on scales from 6 pc djVIueller Sanchez et al] [20061) up to 20 pc 
(IMaiolino et al.lll998l 12000b . This gas most likely is part of the 
cooler regions of the dusty torus. Once clouds enter the ionisa- 
tion cone and are directly irradiated by the central source, they 
emit Ha and [Si vi] emission, which is observed in the ionisation 
cone and which is especially strong on its southern edge. The 
clouds are then accelerated a nd entrained by the outflow (for a 
more detailed discussion, see Packham et al. 2005b- 

As an alternative to the thick toroidal dust distribution in uni- 
fied mode ls, warped disk s are o ften put forward. Fully in line 
with this, IGreenhill et all (120031) have claimed that the warped 
disk alone may be responsible for the obscuration of the nuclear 



3 Note that we are not able to determine an accurate position of the 
MIR source, because high resolution astrometry is not possible with 
MIDI. We hence assume that the torus emission is centered on the 
maser disk, although it might well be possible that there is an offset 
for a slightly inclined torus. 
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source and for the collimation of the outflow in Circinus, obvi- 
ating the need for a geometrically thick dust distribution. Our 
observations, however, have added the decisive component: they 
have proven that a geometrically thick dust distribution is indeed 
present in the nucleus of Circinus. It is obvious that the opening 
angle observed at larger scales is too narrow to be produced by 
merely a disk with a warp of 27°. Hence the disk component can- 
not be the source of the collimation for the outflow and ionising 
radiation alone. 

From our observations we cannot say anything about cooler 
dust, with temperatures T < 100 K, and it is likely that the struc- 
ture we observe is embedded in a larger component of cool gas 
and dust extending out to the starburst on scales of 10 pc. 

6. Conclusions 

We have obtained extensive interferometric observations of the 
nucleus of the Circinus galaxy covering a wavelength range 
from 8.0 to 13.0yum, using the MIDI interferometer at the VLT. 
Through direct analysis of the data and several steps of mod- 
elling with increasing complexity, we have arrived at the follow- 
ing conclusions for the nuclear dust distribution in this active 
nucleus: 

1 . The distribution of dust cannot be described by a single, sim- 
ple component. However, it appears to be distributed in two 
components: 

(a) a dense and warm {T\ > 330 K) disk component with a 
radius of 0.2 pc and 

(b) a less dense and slightly cooler (T2 < 300 K) geomet- 
rically thick torus-like component extending out to 1 pc 
half-light radius. 

2. This disk-torus configuration is oriented perpendicularly to 
the ionisation cone and to the outflow. 

3. The compact, disk-like dust component coincides in orienta- 
tion and extent with the nuclear maser disk. 

4. From the total energy needed to heat the dust and from the 
solid angle intercepted by the dust torus, we infer a lumi- 
nosity of the accretion disk of L acc « 1 • 10 10 L o . This 
corresponds to ~ 20 % of the Eddington luminosity of the 
M BH = 1.7 • 1O 6 M nuclear black hole. 

5. The properties of the larger torus component are less con- 
strained by the data. With an h/r ~ 1, it is consistent with 
a thick, edge-on torus which has cavities for the ionisation 
cone and outflow. This component is most likely responsible 
for collimating the radiation and the wide angle outflow orig- 
inating in the nucleus down to the observed opening angle of 
the ionisation cone and outflow of 90° at larger distances. 

6. We find strong evidence for a clumpy or filamentary dust 
distribution in the torus from three lines of reasoning: 

(a) The large component only has a low effective scaling fac- 
tor of -20%. 

(b) The visibility measurements show an irregular behaviour 
that can be explained by dumpiness. 

(c) The radial temperature dependency of the dust, T(r), has 
a rather shallow decrease which requires that some of the 
outer dust is directly exposed to the nuclear radiation. 

7. The silicate absorption depth is less pronounced for the disk 
component than for the extended component, arguing for 
a 10//m silicate feature in emission in the innermost, disk- 
shaped regions of the dust distribution. 

8. The silicate absorption profile is consistent with that of nor- 
mal galactic dust. We hence see no evidence for dust repro- 



cessing. This might be connected to the moderate maximum 
dust temperatures, T m . dx < 400 K. 

Taking this evidence together, we conclude, that the long pos- 
tulated dust torus exists in this Seyfert 2 nucleus. It is, how- 
ever, significantly more fine structured than the doughnut-shaped 
models which have been proposed as long as no spatially re- 
solved information was available. 
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Appendix A: 1D model 

With this first step of modelling, we aim to estimate the characteristic size of the emission region. To this end, we assumed the 
source to be axisymmetric on the plane of the sky and plotted all visibility measurements, regardless of the position angle, as a 
function of the baseline length BL: Va{BL). This can be done independently for all wavelength bins measured by MIDI. Examples 
for such a plot are shown for two wavelength bins at 8.5 fim and 12.5 fim in the top row of Fig. IA.1I The visibility estimates show a 
decrease towards longer baselines from V = 1 (BL = m) and an irregularly oscillating behaviour for BL > 35 m. At first glance, 
this visibility function is reminiscent of that for a uniform disk, which is a Bessel function of the first kind of order 1: \J\(BL)\. The 
direct comparison in Fig. lA.ll reveals. however, significant discrepancies between such a Bessel function and the data. Additionally, 
our data shows no evidence for the 2n phase jumps that occur at the zero crossings of the Bessel function. We thus conclude that a 
uniform disk model cannot explain the data. 

We chose an alternative approach to analyse the data: First we interpolated Va(BL) (blue curves in the top row of Fig. IA.lt 
and then applied an inverse Fourier transform to this visibility function for each wavelength bin: IJjJ) - Tel [Va(BL)] (d). This 
reconstruction method gives a first order estimate of the radial brightness distribution of the nucleus I,\(d), where d is the angular 
distance from the centre. The shape of the visibility function Va(BL) is only well determined for 20 m < BL < 90 m. This in 
turn implies that we have probed angular sizes d between 10 and 70 mas, considering the entire wavelength range of the MIDI 
observations. These characteristic sizes in the Fourier and the real domain are indicated by vertical dashed and dotted lines in 
Fig. EE 

The reconstructed brightness distribution is composed of two components plus a small fraction of extended flux at distances of 
more than 100 mas of the centre. The extended component is fully resolved by our interferometric set-up, while our measurements 
probe the properties of the other two components. One should bear in mind that we assumed axisymmetry for this analysis and 
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Fig. A.l. Radial ID analysis of the mid-infrared emission in Circinus. Top row: Visibility as a function of baseline length for two 
wavelength bins at 8.5 fim (a) and 12.5 fim (b) taking together all position angles. The data cannot be explained by a uniform disk 
model (light grey curves). The blue curve is an interpolated visibility function derived from the data. Bottom row, (c) and (d): For 
the same wavelengths, the inverse Fourier transform of the blue curve in the top row graphs is plotted, corresponding to a brightness 
distribution. Traced in green is the Gaussian distribution fitted to the Fourier transform in order to get an estimate for the size of the 
emitting source. The vertical lines indicate the angular size range probed by the interferometric measurements. 
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Fig. A.2. Characteristic size of the central dust emission in Circinus as a function of wavelength: the figure shows the FWHM of 
the Gaussian distribution fitted to the radial profile of the brightness distribution of the Circinus nucleus (see Fig. IA. lb . This first 
estimate for the effective size of the nuclear dust emission shows a slow increase with wavelength from 50 to 70 mas. The area 
affected by atmospheric absorption are hatched in grey. 

already from the qualitative discussion of the data in Sect. [3] we know that this is a crude simplification. We therefore used the 
reconstructed brightness distribution only to extract the effective size of the emission region in the mid-infrared and traced its 
evolution depending on the wavelength. This was achieved by fitting a Gaussian distribution to the reconstructed radial profiles at 
every wavelength bin. This Gaussian is shown in green in the bottom row of Fig. IA.ll The full width at half maximum (FWHM) of 
this Gaussian distribution is shown as a function of wavelength in Fig. IA.2l We attribute the small dip in the size between 9.5 and 
10.0 //m to artefacts coming from the ozone feature (hatched area in Fig. |A.2t . From shorter to longer wavelengths, the size of the 
flux distribution increases slightly from 50 to 75 mas. This indicates that the emission at longer wavelengths is more extended than 
at short wavelengths: the temperature of the dust decreases with distance to the nucleus, but apparently only slowly. 

These effective size estimates again demonstrate the need for interferometric observations, as they are the only way to reach 
these resolutions in the infrared. 

Appendix B: 2D geometrical model 

This second model is a two-dimensional, purely geometrical model, aimed to reconstruct the brightness distribution of the source 
from the full two-dimensional visibility data set. The model consists of two concentric, elliptical Gaussian distributions. It was 
motivated by evidence for two components in the one-dimensional analysis and by the dependence of the visibilities on the posi- 
tion angle which imply an elongated source (Sect. |3~2l and there especially Fig. [3}. The brightness distribution has the following 
functional formal 



Here, a and 6 are the position on the sky with respect to the centre of the distribution. This two component model has a total 
of six free parameters: the sizes of the two Gaussians, given by their full width at half maximum Ai{A) and A2(/l); their oblateness, 
represented by the axis ratio ri(A) and riiA); the position angle (p(A) of the major axis with respect to north assumed to be the same 
for both components; and the relative intensity of the two components quantified f(A). 

These parameters were then optimised using the following prescription: for a given wavelength A the flux distribution h{a, $) = 
I(A, a, 6) was generated according to the set of parameters. Then the Fourier transform 'V x(u, v) = T a ,t> [Lxifl, £>)] (u, v) was calculated 
using a fast Fourier transform (FFT). The result corresponds to the visibility distribution in the uv plane. The visibility values for the 
model were extracted from this plane at the positions m, = BL, ■ sin PA,, V; = BL, ■ cos PA, , using the baseline information from our 
observations i consisting of the baseline length BL, and the position angle PA, . Finally, the modelled visibilities V™ od (A) = \V '\{uu v,)| 
were compared to the observed ones V° hs (A). The discrepancy was weighted with the measurement errors cr vl) b S (/V) to obtain the^ 2 : 



4 To be more descriptive, we use the FWHM, A„ = 2 V2 In 2 cr„, in our definition of the Gaussian distribution instead of the sigma, o"„, in the 
standard definition. 
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Fig. B.l. Comparison of the visibilities between the 2D geometrical model (two elliptical Gaussian distributions) and the measured 
data for two wavelengths: (a) 8.5 and (b) 12.5 fim. The data are plotted with asterisks; the modelled visibilities are the continuous 
lines. The two baselines shown here are plotted in two different colours, UT2 - UT3 in green (grey) and UT3 - UT4 in blue (black). 



This x 1 was minimised using the Levenberg-Marquardt least-squares minimisation algorithm. This process was performed indepen- 
dently for the 51 wavelength bins measured by MIDI between 8 and 13 /mi. A comparison of the model to the data is shown for 
8.5 and 12.5 //m in Fig. IB. II This plot is the same as Fig. [4] except that here the visibilities are compared instead of the correlated 
fluxes. The resulting parameter dependencies are shown in Fig. IB. 21 To have more practical parameters, the factor / was replaced 
by the ratio R of the flux in the compact component to the flux in the extended component. The reduced^ 2 (that is ;f 2 (/l)/Mree with 
Affree = 21-6=15 degrees of freedom), shown in the left panel of Fig. IB. 31 varies strongly for different wavelengths. The fits are 
best in the middle of the wavelength range. Towards the edges of the N band, the fits get significantly worse. This is mainly due to 
the weighting of the discrepancy between model and data using the errors of the data cr° hs (A). Due to the bad signal-to-noise ratio 
for low fluxes, the errors of the visibility are of the same order, or even larger within the silicate absorption feature, the atmospheric 
ozone feature, as well as at the edges of the N band. The larger errors lead to a lower x 2 value. When considering a modified estimate 
for the quality of the fit, 



f V° bs (A) - V™ d (A) 



£1*1 v™ d oi) 



(B.3) 



the discrepancy is of the same order, except at A > 12 /im. The reduced £ 2 , i; 2 (A)/Nf lee , is shown in the right panel of Fig. 1531 

With the high values of the reduced x 2 , the fit is comparatively bad. This is partly due to an overly optimistic estimate of the 
errors of the visibilities. The errors only represent the relative errors within one measurement and do not reflect the systematic 
errors between two measurements at different epochs. More importantly, we see in the bad fit quality evidence for a more complex 
structure of the source. This will be discussed in Sect. 14.21 

Although no clear dependency of the parameters with wavelength is identifiable, a few general trends are present. The smaller 
component has a rather constant size of approximately Ai ~ 20 mas and a rather high axis ratio of r\ = 0.2. The position angle has 
a preferred orientation of PA ~ 40°, which is more or less perpendicular to the ionisation cone. Towards both ends of the spectral 
range sampled (8.0//m < A < 9.0 /im and 12.5 //m < A < 13.0 /mi), the position angle deviates significantly from this value. A more 
regular behaviour of the parameters is found in the range from 1 1.0 to 12.5 /mi. Here also the intensity ratio stabilises at R ~ 0.1 
and the size of the second component grows linearly with increasing wavelength. The large component hence comprises the bulk 
of the emission. The ellipticity of the larger component is relatively small, with ri ~ 0.75 for large parts of the spectral range, 
again exceeding these margins at the edges of the spectrum sampled. Some of the parameters show a certain reciprocal influence: 
decreasing the inclination of the smaller component between 8.2 and 9.5 //m leads to a larger emitting region and hence to a larger 
flux, which has an influence on the flux ratio between the two components and the size of the larger component. We hence do not see 
the decrease of ellipticity of the small component to be significant but rather as an artifact of the ambiguity of the data. Especially 
the ellipticity, the size and the flux ratio of the two components show a certain degree of degeneracy. We found that, by fixing some 
of the parameters to average values (e.g. PA to 40° or r\ to 0.2), the quality of the fit for the individual wavelengths does not decrease 
significantly. At the same time, however, the behaviour of the other parameters, which are then adjusted, looks steadier. We did not 
follow this path of analysis any further but decided to apply a physical model instead (see Sect. 14.2b . The overall erratic behaviour of 
the parameters and the quality of the fit, however, show that two Gaussians are not a good fit and that there are substantial differences 
at different wavelength bins. Nonetheless, we find that other models of the same simplicity produce no better fits. A more detailed 
discussion of alternative models is be given in Sect. 14.11 
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Fig. B.2. Wavelength dependence of the parameter estimates for the two-dimensional geometric model. The six parameters are the 
two FWHM of the Gaussian components (a), the respective axis ratios r (b), the position angle <p of the elliptical distribution (c) and 
the ratio R of the flux in the compact over the extended component (d). The areas affected by atmospheric absorption are hatched in 
grey. 
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Appendix C: Comparison of the full data set to the physical model 
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Fig. C.l. Comparison of the dispersed correlated fluxes (green/gray) to the model fluxes for the physical model fit (fit2, blue/black) 
in the entire wavelength range from 8 to 13 /an. The dashed line represents the flux contribution of the small disk component, while 
the dotted line is the contribution of the large torus component. 



